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Recent experiments show that indirect resonant inelastic X-ray scattering (RIXS) is a new probe 
of spin dynamics. Here I derive the cross-section for magnetic RIXS and determine the momentum 
dependent four-spin correlation function that it measures. These results show that this technique 
offers information on spin dynamics that is complementary to e.g. neutron scattering. The RIXS 
spectrum of Heisenberg antiferromagnets is calculated. It turns out that only scattering processes 
that involve at least two magnons are allowed. Other selection rules imply that the scattering 
intensity vanishes for specific transferred momenta q, in particular for q = 0. The calculated 
spectra agree very well with the experimental data. 



Introduction. Resonant inelastic X-ray scattering 
(RIXS) is a technique that is rapidly developing due to 
the recent increase in brilliance of the new generation 
synchrotron X-ray sources [H . As a scattering technique 
RIXS has two important advantages. First, it is sensitive 
to excitations that are difficult to observe by otherwise, 
for example direct d-d excitations in cuprates or mangan- 

ites [l[a0,[!,@,0]- Second ' it: probes directly both the 
energy and momentum dependence of such excitations 
-unlike for instance optical methods. 

Experiments are performed at the K-edges of transi- 
tion metal ions, where the frequency of the incident hard 
X-rays is tuned to match the energy of an atomic ls-Ap 
transition (5-10 keV). When a photon is absorbed in the 
solid a Is core- hole is created as the electron is promoted 
to the 4p empty states. From the theoretical side it is 
clear that the scattering intensity is related to the charge 
dynamics of the system under study [H, 0, [l(J, 11 | ■ The 
precise nature of this relation was established recently: 
a systematic expansion in the ultra-short life time of the 
core-hole [12J, |l3| shows that the effective RIXS cross sec- 
tion a resonant scattering factor times a combination of 
the linear charge response function S(q,uj) and the dy- 
namic longitudinal spin density correlation function of 
the d— electrons. 

A recent breakthrough was achieved by J. P. Hill et 
al. 14j, who observed that RIXS on the high temper- 
ature superconductors L^-zCu^C^ picks up transver- 
sal spin dynamics -magnons. This triggers the question 
what kind of transversal spin correlation function actu- 
ally is measured in magnetic RIXS and how this is dif- 
ferent from other probes of spin dynamics -in particular 
inelastic neutron scattering. The aim of this Letter is to 
answer that question. 

Experimentally the magnetic RIXS signal is already 
already in the undoped cuprate La2Cu04 -in fact it is 
strongest in this case. The magnetic loss features are at 
energies well below the charge gap of this magnetic in- 
sulator, at energies where the charge response function 
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FIG. 1: Schematic representation of the magnetic RIXS scat- 
tering process at a transition metal K-edge. Left: the incom- 
ing photon (energy wf n , momentum qi n ) induces an electronic 
transition from a Is to a 4p level. Middle: exchange inter- 
action between 3d electrons in the presence of the core-hole. 
Right: de-excitation and outgoing photon (Wout 5 Qout)- 



S(q,uj) vanishes, as well as the longitudinal spin one - 
which is in fact a higher order charge response function. 
This implies that we need take the ultra-short lifetime 
12, [lj| one step further. By doing so I find 



expansion 

the scattering amplitude for magnetic RIXS, expressed 
in terms of an intrinsic dynamic four-spin correlation 
function of the system that is probed. Moreover I de- 
rive selection rules which are related to the symmetry of 
the underlying spin Hamiltonian. It turns out that the 
first allowed magnetic scattering process is a two magnon 
scattering one. As an example the formalism is used to 
calculate the indirect RIXS spectrum of Heisenberg an- 
tiferromagnets as function of transferred momentum q. 
The scattering intensity vanishes at (0, 0) and at the an- 
tiferromagnetic wavevector q = (jr, 7r). The computed 
spectra all agree with the experimental data. 

Series expansion of the scattering cross section. The 
probability for X-rays to be scattered from a solid state 
system can be enhanced by orders of magnitude when 
the energy of the incoming photons is in the vicinity of 
an electronic eigenmode of the system -i.e. in the vicin- 
ity of a resonance. At a transition metal K-edge a Is 
electron from the inner atomic core is excited into an 
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empty 4p state, see Fig. [Tj In transition metal systems 
the empty 4p states are far (10-20 eV) above the Fermi- 
level, so that the X-rays do not cause direct transitions of 
the Is electron into the lowest 3<i-like conduction bands 
of the system. Still this technique is sensitive to low en- 
ergy excitations of the d— electrons because the Coulomb 
potential of Is core-hole can couple to e.g. very low en- 
ergy electron-hole excitations when the system is metal- 
lic. Since the charge excitations are caused by the core- 
hole, this scattering mechanism is sometimes referred to 
as indirect RIXS. 

In this Letter, however, we will consider insulating sys- 
tems -in particular Mott insulators, where the only re- 
maining low-energy degrees of freedom are the spin ones. 
So in order to determine the RIXS scattering amplitude 
we need to establish how in this case the core-hole cou- 
ples to magnetic excitations. Our starting point is the 
Kramers-Heisenberg formula for the resonant scattering 
cross section 15, 16] 
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and / and i denote the final and initial state of the 
system, respectively. The sum is over / is over all fi- 
nal states. The momentum and energy of the incom- 
ing/outgoing photons is qj n /out and Wj° n y out and the loss 
energy u> — uj® ut — u>f n is equal to the energy difference 
between the final and initial state u>fi — Et — E^. In 
the following we will take the groundstate energy of our 
system as reference energy: Ei = 0. In the scattering 
amplitude Afi the resonant energy is w rcs , n denotes the 
intermediate states and D is the dipole operator that de- 
scribes the excitation from initial to intermediate state 
and the de-excitation from intermediate to final state. 
The dipole operator is given in more detail in for in- 
stance Ref . 12| . The energy of the incoming X-rays with 
respect to the resonant energy is u>i a (this energy can thus 
either be negative or positive: ui- m = uj® n —ui Ies ) and E n is 
the energy of intermediate state \n) with respect to the 
resonance energy. The last important detail is that the 
intermediate state is not a steady state. The highly ener- 
getic Is core-hole quickly decays e.g. via Auger processes 
and the core- hole life-time is very short. This leads to a 
core-hole energy broadening Y which is proportional to 
the inverse core-hole life-time. 

To calculate RIXS amplitudes, we proceed by formally 
expanding the scattering amplitude in a power series [121 ] 
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where we introduced A = uj- ln - 
in the intermediate state H lrA . 



iT and the Hamiltonian 
For a further expansion 



of this scattering amplitude it is essential that we split 
up the intermediate state Hamiltonian into two parts: 
Hi n t = Hq + Hi, where Ho is the Hamiltonian of the 
system without core-hole and Hi the part of the Hamil- 
tonian that is active in the presence of a core-hole. 

Spin Hamiltonian with core-hole. We will calculate the 
resonant X-ray cross section in a Mott-Hubbard insulator 
at zero temperature. We assume that this system is de- 
scribed a single band Hubbard model at strong coupling 
and at half filling. In this case the electrons are local- 
ized and the only low energy degree of freedom is their 
spin. It is well known that in a Mott-Hubbard insula- 
tor the magnetic exchange integrals are determined by a 
virtual hopping process of electrons. We denote the hop- 
ping amplitudes of the valence electrons by ty where i 
and j denote lattice sites with lattice vectors R; and Rj . 
The Coulomb interaction between electrons at the same 
site is U, so that in second order perturbation theory we 
have the exchange interaction Jjj = 2tfj/U and the spin 
dynamics is governed by a Heiscnbcrg spin Hamiltonian 
of the form 
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where Jk is the Fourier transform of . It is well known 
that for neutron scattering on such a spin system the two- 
spin correlation function J2 a I e~ iuJt (Sq (Q)S-q(t))dt is 
measured, where the sum a is over the three spin compo- 
nents. We will see shortly that magnetic RIXS measures 
a very different four-spin correlation function. 

In the intermediate state a core- hole is present. We as- 
sume the core-hole potential U c to be local, i.e. as acting 
exclusively on those (valence) electrons that belong to 
the atom with the core-hole. When on site m a core-hole 
is present, the exchange interactions that involve the spin 
on site m becomes stronger, as the virtual intermediate 
state with two electrons on the site with the core-hole 
are lowered in energy by U c . On the other hand the vir- 
tual state with two holes present on site m is at U + U c . 
Adding these two effects leads to the following exchange 
interaction between the spins on site m and j of 
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and 77 = jjtzjJ2- From this we obtain Hi, part of the 
intermediate state Hamiltonian that is active in the pres- 
ence of a core-hole 
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where the operator s m creates a core-hole on site m. 

Spin-spin correlation function as measured in RIXS. 
In order to finally obtain the magnetic cross section, we 
need to evaluate the operator (-ffint)' in equation 
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FIG. 2: Left: RIXS spectrum for a nearest neighbor Heisen- 
berg antiferromagnet with exchange interaction J as a func- 
tion of transferred momentum q for a cut through the 
Brioullin zone that is show on the right. 



This is a non-trivial task. We first expand (i?i n t) m 
a series that contains the leading terms to the scat- 
tering cross section in lowest order in rjJ/ A. A con- 
servative estimate gives, using for the copper K-edge 
r Ri 1.5 eV and for the high temperature superconduc- 
tors J w 125 meV and U c /U « 0.85 @,[Lz|], that at reso- 
nance r]J/A w 0.22. This makes it a suitable expansion 
parameter. Re-summing the leading order terms in the 
series gives in the end the magnetic scattering amplitude 
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where we find that the scattering operator O q to be 
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Hamiltonian of the system (O q =o ^ Hq). Thus |z) and 
\f) are eigenstates of this magnetic scattering operator, 
which makes inelastic scattering impossible [18j. This is 
in stark contrast with conventional two-magnon Raman 
scattering in the optical or UV range. That technique 
is also sensitive to a four-spin correlation function [19(, 
but a quite different one. This is obvious considering the 
fact conventional Raman scattering is restricted to q = 
-precisely the momentum transfer where RIXS vanishes. 
Therefore also these two techniques offer complementary 
information on spin dynamics. 

Two-magnon scattering in antiferromagnets. From 
equation ([7]) immediately another selection rule follows. 
The projection of the total spin on the z-axis, S^ ot = 
J2 i Sf commutes with both H and O q . Therefore Sf ot 
is conserved during the scattering process, which implies 
that creation of a single magnon by the core-hole is not 
possible. But the creation of two magnons (with opposite 
z-projections) is allowed and this is therefore the lowest 
order transversal spin scattering process that contributes 
to indirect RIXS. Also four-magnon scattering is in prin- 
ciple allowed, but of higher order and therefore smaller 
and not taken into account in the linear spinwave analy- 
sis that follows. Note that on the grounds of symmetry 
it is possible, in principle, to have magnetic scattering 
without creating any additional magnons in the scatter- 
ing process. Physically this situation can only arise at 
finite temperature, when a magnon with momentum k 
that is present in the groundstate is scattered to k + q 
by the core-hole. This implies that magnetic RIXS has an 
interesting temperature dependence -which is, however, 
beyond our present scope. 

We now apply the theory above to two-dimensional 
bipartite S = 1/2 antiferromagnets and determine the 
two-magnon RIXS spectrum as a function of transferred 
momentum, at zero temperature. To this end the Hamil- 
tonian Hq and correlation function O q are bosonized 
within linear spinwave theory, where — > a~l , — » as, 
and Sf — > \ — n.i, with boson creation/annihilation op- 
erators ai/df and number operator m = afa,i. After 
a Boguliobov transformation into the boson operators 
cti/af we have a J = u^a^ + Ukfl-k with 



so that the magnetic correlation function that is mea- 
sured in RIXS is proportional to the spin correlator 
Je-^<<5 q (0)d_q(t)>dt. 

This expression is surprisingly simple and elegant. It 
shows that indeed momentum resolved indirect RIXS 
probes a momentum dependent four-spin correlation 
function. From expression it is immediately clear why 
experimentally the magnetic RIXS intensity vanishes at 
zero transferred momentum, i.e. at q = 0. In that case 
the correlation function is nothing but the steady state 



and ek = 2 y/ J^ =Q — , then the Hamiltonian reduces to 
jjLSW _ ^ k e k a^a:k. It is straightforward to show now 
that within linear spinwave theory the two magnon part 
of the magnetic scattering operator is 
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FIG. 3: Top: first moment of the RIXS spectrum and bot- 
tom: total spectral weight for a nearest neighbor Heisenberg 
antiferromagnet (full line) and a Heisenberg model with in 
addition a second and third neighbor exchange (J' and J", 
respectively) (dashed line). 



where q is the total momentum of the two magnon excita- 
tion. The resulting RIXS spectrum is shown in Fig. [2j for 
a cut through the Brioullin zone indicated by the right- 
hand side of the figure. There are several remarkable 
features in the spectrum. 

First of all the spectral weight vanishes at q — (0, 0) 
and q = (tt,tt). This in agreement with the experimen- 



tal observations [14| ■ From the scattering operator ([9]) it 



one sees that this selection rule is due to the antiferro- 
magnetic ordering that occurs in the Heisenberg Hamil- 
tonian. It is actually easy to show that the RIXS inten- 
sity always vanishes at (7r, 7t) if this scattering vector is 
also a reciprocal lattice vector. This holds for instance 
also for a Heisenberg Hamiltonian with weak second and 
third neighbor exchange intereactions (J' and J", respec- 
tively), which is illustrated by the calculation RIXS spec- 
trum for an extended Heisenberg antiferromagnet with 
J' = J" = J/20, shown in Fig. [3l The longer range 
couplings transfer spectral weight to scattering vectors 
around (ir, it), but do not induce weight at precisely that 
wavevector. 

The other remarkable feature feature of the magnetic 
RIXS spectrum is its strong dispersion. This is apparent 
from Fig. [5] and the upper pannel of Fig. O which shows 
the first moment (average peak position) of the spectrum. 
The calculations for the nearest neighbor Heisenberg an- 
tiferromagnet (Fig. [31) show that the magnetic scattering 
disperses from about oj k, around (0, 0) and (tt, tt) to 
w w 4 J at (w, 0) and (tt/2, 7r/2). Longer range couplings 
tend to reduce the first moment of the RIXS spectrum. 
The observed dispersion has a two-fold origin. It is in 
part due to the q-dependence of the two-magnon density 
of states (DOS), combined with the scattering matrix el- 
ements that tend to pronounce the low energy tails of the 
two-magnon DOS. 

The consistency at q = (0, 0) and q = (tt, tt) of the 
theoretical results and experimental data was already 



noted, but at other wave-vectors the agreement stands 
out even more. The data on La 2 Cu04 show for q = (tt, 0) 
a peak at around 500 meV, precisely where we find it on 
the basis of a nearest neighbor Heisenberg model with 
J = 146 meV -a value also found by the analysis of neu- 
tron scattering data [2(J. Similar agreement is found at 
q = (0.6tt, 0) and q = (0.6tt, 0.6vr). 

Conclusions. We determined the momentum depen- 
dent four-spin correlation that is measured in magnetic 
RIXS. On the basis of this the magnetic RIXS spec- 
trum was calculated for Heisenberg antiferromagnets 
with short and longer range couplings. We derive selec- 
tion rules show that only scattering processes that involve 
at least two magnons are possible and that the scatter- 
ing intensity vanishes at zero momentum transfer and at 
the antiferromagnetic lattice vector -which are observed 
in experiment. Moreover theory and experiment agree 
very well on the two-magnon peak position and spectral 
weight throughout the measured Brioullin zone. These 
results show that RIXS is in principle a powerfull tool to 
obtain new information on spin dynamics -information 
that is complementary to what can be obtained by other 
techniques such as neutron, non-resonant X-ray or con- 
ventional two-magnon Raman scattering. 
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